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Introduction 

Due  to  very  large  curvature  of  both  sides  of  a  typical  turbine  blade,  the  flow  past 
each  of  them  has  specific  features  not  intrinsic  to  wing  aerodynamics.  The  Figure  below 
illustrates  this  statement,  showing  a  long  region  of  laminar/turbulent  transition  on  the 
upper  (suction)  side  and  absolutely  unstable  Gortler  vortices  adjacent  to  the  lower 
(pressure)  side.  The  streamwise  pressure  gradient  is  favourable  over  most  of  the  suction 
side  and  adverse  on  the  pressure  side  that  provokes  turbulent  separation.  Thus,  the  types  of 
instabilities  to  be  expected  to  develop  on  each  side  should  be  different  in  nature.  Turbulent 
separation  under  the  influence  of  the  strong  adverse  pressure  gradient  is  an  additional 
poorly  understood  phenomenon  characteristic  of  turbomachinery  flows.  Up  to  recently,  all 
these  topics  were  beyond  the  scope  of  any  rigorous  mathematical  treatment,  challenging 
theoreticians  to  develop  a  pertinent  high-Reynolds  number  asymptotic  approach  which 
would  cover,  along  with  instabilities,  also  late  (or  deep)  transition.  Traditionally,  transition 
is  regarded  as  the  most  intricate  stage  of  flow  on  the  turbine/compressor  blade  even  as 
compared  to  fully  developed  turbulence.  However,  turbulent  velocity  fields  are  hardly 
amenable  to  analytical  study.  Only  a  few  attempts  were  made  to  attack  the  simplest  case  of 
the  zero-pressure-gradient  flow  on  a  flat  plate  without  introducing  special  assumptions  on 
the  closure  model.  Previous  theoretical  studies  of  turbulent  separation  involved  crude 
errors.  Thus,  the  proposed  research  is  aimed  at  resolving  long-standing  difficult  problems 
of  fundamental  importance  from  both  scientific  and  industrial  standpoints,  especially  as 
applied  to  gas-turbine  engine  aerodynamics. 

More  specifically,  the  work  done  under  the  grant  started  on  January  1,  2003 
encompasses  the  following  topics: 

1 .  Gortler  vortices  on  the  concave  pressure  side  of  a  blade; 

2.  Laminar/turbulent  transition  on  the  convex  side  of  a  blade; 

3.  Turbulent  separation  under  the  action  of  the  strong  adverse  pressure  gradient; 

4.  Cooling-film  aerodynamics. 

In  fact,  the  cooling-film  flows  fall  beyond  the  scope  of  the  proposal,  but  a  shift  in  priorities 
in  transition  on  gas-turbine  blades  made  it  timely  to  extend  the  contents  of  the  proposed 
research  so  as  to  cover  this  new  important  issue.  Note  that  to  simplify  notation,  a  separate 
numbering  of  equations  and  figures,  as  well  as  a  separate  list  of  references,  is  used  in  each 
of  the  sections  below. 
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Suction  Side 


A  long  region  of 
laminar/turbulent  transition 
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j  /  \\rrheory:  Ryzhov  2005) 


Pressure  Side 

Absolutely  unstable  Gortler  vortices 
{Ryzhov  &  Bogdanova-Ryzhova  2004) 

Figure. 


Part  I:  Gortler  vortices 


As  it  has  been  proposed,  the  formulation  with  the  vibrating  ribbon  was  chosen  for 
simulating  basic  properties  occurring  in  real  turbomachinery  environment  as  well  as 
experimental  setups. 


1.  Receptivity 


Within  the  framework  of  an  asymptotic  approach  based  on  the  triple-deck  theory, 
the  initial  system  of  Navier-Stokes  equations  reduces  in  the  viscous  near- wall  sublayer  to  a 
much  simpler  set  of  Prandtl  equations  (  Stewartson  1974  ) 
du  +  dv  dw  _ 
dx  dy  dz 


du  du  du 

—  +  u +  V —  + 

dt  dx  dy 


du  dp  d2u 

w —  = - +  — t- 

dz  dx  dy 


(la,b,c) 


dw  dw  dw  dw  dp  d2w 

- +  U +  V —  +  w — •  =  -  —  + — 

dt  dx  dy  dz  dz  dy 


for  an  incompressible  boundary  layer,  no  matter  what  the  Mach  number  of  the  oncoming 
flow.  The  pressure  p=  p{t,x,z)  not  known  in  advance  from  a  solution  of  the  Euler 
equations  is  to  be  determined  simultaneously  with  the  velocity  field  and  the  instantaneous 
displacement  thickness  -A(t,x,z) .  Crossflow  vortical  disturbances  obey 
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whereas  unsteady  Gortler  vortices  are  controlled  by  (  Ryzhov  &  Bogdanova-Ryzhova 
2003) 
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Here  zD  and  s  3H  are  similarity  parameters,  the  wall  curvature  k  is  included  in  H,  z  -»  0 
as  the  Reynolds  number  Re  ->  oo .  The  size  of  disturbances  is  fixed  by  the  pressure 
variations  which  are  AjV2)  while  the  amplitude  of  a  perturbing  agency  is  as  small  as 

djz5 ) .  When  dealing  with  an  orderly  roughness,  we  introduce  the  time-dependence  by 
means  of  a  factor  exp(zco0r)  for  t>  0  and  assume  the  periodic  obstacle  geometry  to  be 
specified  as 

yw (v, z)  =  q/(x) cos(^z),  t>  0  (4) 


In  this  approach  yw(x,z)  =  0  for  /  <  0 .  The  main  focus  is  on  disturbances  emitted  during 


the  initial  pulse  mode  at  t  close  to  0  when  the  generation  of  highly  modulated  wave 
packets  rather  than  harmonic  Tollmien-Schlichting  wave  trains  takes  place.  This  is  the 
crucial  point  of  the  mathematical  model  since  the  most  amplifying  wave  packets  propagate 
with  the  group  velocity  as  distinct  from  the  phase  velocity  determining  the  direction  in 
which  the  Tollmien  -  Schlichting  waves  sweep  downstream  of  a  perturbing  source 
(Landau  &  Lifshitz  1959  ).  The  group  and  phase  velocities  are  different  in  general; 
accordingly,  the  modulated  wave  packets  and  periodic  wave  trains  can  be  driven  in  the 
opposite  direction. 

On  the  assumption  that  a  parameter  a  entering  (4)  tends  to  0,  a  solution  to  the 
receptivity  problem  based  on  the  linearized  equations  (la-c)  reads 
p=  Wj ocem'"\  (4a) 
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with  S  standing  for  either  of  the  two  similarity  parameters  zD  or  z  'H  depending  on 
whether  the  crossflow  or  Gortler  vortices  are  under  consideration.  Here  O(d)  is  a  well- 
known  function  expressed  through  the  first  derivative  d4i{fl)/dY  and  an  improper  integral 

oo 

/(d)  =  J  Ai(  Y)dYoi  the  Airy  function  Ai(  Y)  by  means  of 
n 

0(n)  =  ^j&r'(n)  (5) 

The  argument 
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comprises  the  frequency  co  ,  both  wavenumbers  k,  and  the  components  xx,  x,  of  the 
wall  shear  stress.  For  wave  system  on  a  concave  wall,  we  have  xz  =  0,  k=  1  and  the 
generalized  wavenumber  simplifies  to  K=k.  In  this  case  the  right-hand  side  of  the 
dispersion  relation 


O(Q)  =  £(k,m0) 


(6) 


takes  the  form 
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having  a  strong  singularity  k  3 .  However  the  generalized  wavenumber  K  determines  an 
analogous  singularity  in  the  right-hand  side 
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of  the  dispersion  relation  (6)  for  crossflow  vortices  governed  by  the  interactive  law  (2). 

2.  Residue  series  expansion 

With  real  k,  both  positive  and  negative,  the  complex  frequency  plane  has  no 
branch  points.  Applying  residue  theorem  to  the  Laplace  integral  in  (4b)  yields 
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where  ;Q0  =  P(a0K  3  and,  in 
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keeping  with  above  abbreviations, 
are  the  dispersion  relation  roots.  The  first 
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term  A,  in  (9a)  is  periodic  in  time  and  has  nothing  to  do  with  exponentially  growing  wave 
packets  which  are  of  primary  concern  for  our  purposes.  As  is  known  from  the 
hydrodynamic  stability  theory,  only  the  first  root  co ,  in  (9c)  gives  rise  to  unstable 
oscillations  whereas  all  the  other  eigenmodes  related  to  the  roots  co ,(&),  j  =  2,3,.. ,  damp 

out  with  time  both  upstream  and  downstream  of  the  perturbing  source.  In  an  asymptotic 
representation  for  A,  we  may  confine  ourselves  to  the  leading-order  approximation  and 
write 
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3.  Reduction  to  single  integrals 


With  (10)  in  hand,  an  inverse  Fourier  transform  designated  by  ma, in 
(4b) reduces  to 
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for  crossflow  vortical  disturbances  and 
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for  unsteady  Gortler  vortices.  A  study  of  the  single  integrals  obtained  presents  a 
mathematical  problem  not  nearly  as  complicated  as  the  initial  one  given  in  (4b).  Rigorous 
results  based  on  (1 1)  and  (12)  are  set  forth  below. 


4.  Spectrum 

In  the  Tollmien-Schlichting  interval  of  frequencies  and  wavenumbers  the  first 
terms  in  curled  brackets  on  the  right-hand  sides  of  (7)  and  (8)  are  dominant.  The  correction 
terms  become  comparable  in  magnitude  to  the  leading-order  terms  within  two  spectral  side 
bands.  A  high-frequency  side  band 
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"k  =  (s/T7  mcrfi  ~  2-0 mcrfi 

specifies  the  range  of  unsteady  crossflow  vortices  (  estimates  here  relate  to  a  transitional 
boundary  layer  with  Reynolds  number  Re  ~5  •  1 0s ).  A  low-frequency  side-band 

\_ 

CO  =(s6//-)7cofJ  ~0.25©(; 

kG*0.\2kG  (14) 

”h  =  [^H  )  7  m<;,o  ~  2.0/t iGfi  «  2.0wtr  0 

defines  the  range  of  unsteady  Gortler  vortices  (  estimates  are  made  with  the  same 
transitional  Reynolds  number  ).  Notice  that  the  vortices  of  both  types  are  equally  packed  in 
the  spanwise  direction  and  may  be  hardly  discernable  from  each  other  if  the  three- 
dimensional  (swept)  and  surface  curvature  effects  simultaneously  play  an  important  role  in 
creating  centrifugal  forces  to  balance  out  the  normal-to-wall  pressure  gradient.  Thus,  the 
disturbance  spectrum  consists  of  the  two  side  bands  with  the  Tollmien-Schlichting  interval 
of  eigen  frequencies  in  between. 

5.  Phase  and  group  velocities 


It  is  obvious  from  the  large-time  solution  to  the  receptivity  problem  posed  in  (4) 
that  the  overall  system  of  pulsations  is  determined  by  the  behavior  of  the  first  root  co ,  of 
the  dispersion  relation  (6)  with  the  right-hand  side  prescribed  either  by  (7)  or  by  (8).  The 
phase  velocity 

c=~  (15) 

ik 


of  the  monochromatic  Tollmien-Schlichting  wave  train  differs  from  the  group  velocity 

i/  _  gh)  , 


evaluated  at  the  points  where  3f(co,)  attains  its  global  and  local  maxima.  Since  wave 


packets  propagate  at  the  group  velocity,  nothing  can  be  said  about  the  direction  of  their 
motion  from  (15)  as  it  is  sometimes  done  when  predicting  transition. 


6 


Figure  1 


Figure  1  exhibits  the  complex  frequency  plane  for  Tollmien-Schlichting  waves 
(using  the  Gortler  variables  with  mG0  =  0).  The  first  dispersion  curve  coG1 ,  responsible  for 
inducing  unstable  disturbances,  consists  of  two  branches  positioned  symmetrically  about 
the  real  axis.  Each  branch  involves  the  only  lobe  with 

^Kl)<0  (17) 

dk(i 

at  every  point.  In  this  simplest  case  the  directions  of  the  phase  and  group  velocities 
coincide  with  each  other.  The  Blasius  boundary  layer  offers  an  example  marked  by  the 
conventional  scenario  of  convective  instability. 

6.  Mechanisms  underlying  streamwise  absolute  instability 


The  key  issue  is  what  happens  to  the  boundary-layer  stability  properties  under  the 
action  of  centrifugal  forces.  The  steepest  descent  method  is  best  suited  to  solving  this 
problem.  The  existence  of  saddle  points  of  the  phase  function 
<p  =  w0)  +  iVsk,  Vs  =  x/t  entering  the  exponent  of  the  integrand  in  (11)  and  (12)  is 

central  to  this  method;  their  coordinates  ks  =  k^n^,  come  from 


dk 


+  iVs  =  0 


(18) 


A  complete  set  of  saddle  points  is  needed  to  compute  the  pressure  variations  along  an 
arbitrary  ray  x/t  =  const .  It  should  be  recognized  that  a  saddle  point  k-ks  contributes  to 
instability  if  the  original  integration  path  over  the  real  axis  can  be  continuously  deformed 
into  a  steepest  descent  contour  through  this  particular  point.  However  we  may  put  Vs  =  0  in 
(18)  to  settle  the  most  significant  issue  of  absolute  instability  in  the  streamwise  direction. 
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The  saddle  point  of  the  phase  function  specified  by  this  condition  becomes  simultaneously 
a  saddle  point  of  the  complex  frequency.  In  numerous  checks  executed  the  original 
integration  path  in  the  Fourier  transforms  (11)  and  (12)  was  found  to  be  deformable  into 
the  steepest  descent  contours  through,  respectively,  the  saddle  points  cocH  =  cotrl(/cu.,  wtr0) 


and  <Bf;i  =  coGi(^;, w(;0) .  Thus  we  are  left  to  conclude  that  absolute  instability  in  the 

streamwise  direction  of  the  boundary  layer  is  triggered  by  the  vortical  (crossflow  or 
Gortler )  disturbances  in  response  to  the  initial  pulse-like  motion  of  the  ribbon. 

A  close  inspection  of  the  shapes  of  the  first  dispersion  curve  in  the  crossflow  and 
Gortler  complex  frequency  planes  reveals  the  mechanism  of  the  streamwise  absolute 
instability.  Figure  2  shows  that  the  first  dispersion  curve  of  crossflow  oscillations  falls 
apart  into  two  branches  located  asymmetrically  about  the  real  axis.  Each  branch  consists  of 
two  lobes.  Both  right-hand  lobes  having 


fgKJ 

dkcr 
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(19) 


are  at  the  bottom  of  the  conventional  scenario  of  convective  instability  whereas  the  two 
left-hand  lobes  form  a  new  crossflow  eigenmode  with 
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at  any  point. 
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Figure  2 


This  vortical  eigenmode  controls  oscillations  in  front  of  a  site  where  they  are  given 
birth.  What  is  more,  the  new  eigenmode  is  endowed  with  a  local  peak  do f  iH(cotr)  at  the 
tip  of  the  loop  that  connects  the  right-hand  and  the  left-hand  lobes  of  the  branch  situated  in 
the  lower  half-plane  of  the  complex  coc, .-plane.  This  local  peak  is  responsible  for  emitting 
modulated  wave-packets  upstream,  against  the  oncoming  flow  ( Landau  &  Lifshitz  1959). 
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Figure  3 


Figure  3  drawn  with  kG  >0  on  an  enlarged  scale  illustrates  the  formation  of  the 
loop  in  the  lower  branch  of  the  first  dispersion  curve  for  Gortler  vortical  oscillations  in  the 
high  subsonic  Mach  number  regime.  The  upper  branch  corresponding  to  kG  <  0  is  a  mirror 
image  of  the  lower  branch.  As  is  obvious  from  this  plot,  the  phase  velocity  c(&)  defined  by 

(15)  is  positive  independent  of  the  sign  of  k.  Hence  the  harmonic  wave  trains  propagate 
always  downstream  triggering  convective  instability  of  the  boundary  layer  on  a  concave 
surface.  On  the  contrary,  modulated  signals  driven  by  the  local  peak  <7  of  93(co  f;)  and  a 
similar  peak  on  the  upper  branch  advance  upstream  with  the  group  velocity  ^<0 

according  to  (16).  Thus  it  may  be  reasonably  predicted  from  the  curves  in  Figures  2  and  3 
that  exponentially  amplifying  oscillations  develop  not  only  in  the  immediate  vicinity  of  a 
perturbing  agency  but  move  upstream  of  it  in  the  form  of  wave  packets  consisting  of  long- 
scaled  cycles  (small  k ).  Upstream  advancing  signals  underlie  the  mechanism  of 
streamwise  absolute  instability  provoked  by  the  wave  /  vortex  eigenmodes’  interaction. 
Mathematically,  the  two-lobe  structure  of  either  of  the  two  branches  of  the  first  dispersion 

_5 

curve  is  brought  about  by  a  singularity  k  3  in  the  right-hand  side  of  the  dispersion  relation 

_5 

(7)  and  an  analogous  singularity  K  3  involved  in  (8). 

7.  Computed  results 

In  line  with  the  above  theoretical  predictions,  the  wave  packets  are  seen  in  Figures 
below  to  form  both  downstream  and  upstream  of  the  vibrating  ribbon.  No  special  thermal 
source  was  included  in  the  computations  where  high  frequency  cycles  were  partially 
filtered  out.  The  high  subsonic  Mach  number  regime  brings  out  the  significance  of  the 
streamwise  absolute  instability  caused  by  the  wave  /  vortex  eigenmodes  interactions  at  any 
value  of  A4,  <  1  • 
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Figure  4 

Figure  4  shows  the  wave  system  for  the  boundary  layer  with  crossflow  specified  by 
mcr0  =  0.4 1 7  and  =  0.95 .  The  first  dispersion  curve  in  the  complex  crossflow  frequency 
plane  is  depicted  in  Figure  2  for  this  case.  Figure  5  gives  the  disturbance  pattern  for  the 
boundary  layer  on  a  concave  surface  that  corresponds  to  the  first  dispersion  curve  in  the 
complex  Gortler  frequency  plane  of  Figure  3.  The  downstream  moving  wave  packets 
underlie  the  conventional  route  to  convective  instability.  The  upstream  advancing  wave 
packets  are  at  the  heart  of  streamwise  absolute  instability  developing  in  front  of  a 
perturbing  source.  As  predicted  theoretically,  streamwise 
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Figure  5 

absolute  instability  is  brought  about  by  the  new  crossflow  and  Gortler-vortex  eigenmodes 
endowed  with  positive  values  of  c£s(<s) , )  jdk . 
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Part  II:  Laminar /turbulent  transition  on  the 
convex  side  of  a  blade 


1  Physical  reasoning 

External  perturbing  sources  operating  in  the  pulse  mode  generate  strongly- 
modulated  wave  packets  in  the  Tollmien-Schlicting  spectral  range  which 
cause  transition  to  complete  long  before  reaching  the  upper  stability  branch. 
The  continuous  filling-up  of  the  wavenumber  spectrum  by  the  overlapping 
wings  of  neighbouring  peaks  gives  rise  to  erratic  short-scaled  wiggles  cor¬ 
rupting  the  primary  long-scaled  oscillation  cycles.  In  this  regard  transi¬ 
tion  induced  by  the  time-harmonic  agency  shows  marked  distinctions  from 
the  explosive  transitioning  to  turbulence  experienced  by  the  wave  packets. 
This  simple  preliminary  consideration  underlies  the  analysis  below  where  the 
triple-deck  formalism  is  exposed  in  more  detail  than  in  section  1.  The  ex¬ 
tension  is  aimed  at  providing  the  clear  understanding  of  the  reasons  why  the 
powerful  asymptotic  approach  that  holds  for  laminar  boundary  layers,  as  it 
was  initially  designed,  can  be  extended  to  transitional  regimes. 


2  Theory 

Reasoning  from  the  experimental  findings  of  Gostelow,  Hodson  &  Walker 
(1999)  and  Gostelow  &  Thomas  (2005)  that  instabilities  develop  in  the 
Tollmien-Schlichting  spectral  range,  we  begin  with  the  triple-deck  scheme 
set  forth  in  Stewartson  (1974).  This  approach  leans  upon  asymptotic  expan¬ 
sions  in  a  small  parameter 


e  =  Re  »  , 


(2.1) 
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with  Re  being  the  characteristic  Reynolds  number.  When  considering  exper¬ 
imental  data,  Re  is  as  a  rule  evaluated  with  the  distance  L*  counted  along 
a  section  from  its  nose.  The  empirical  correlation  by  Mayle  (1991)  gives  the 
transition  length  Reynolds  number,  Re^,  for  the  blade  suction  side  in  terms 
of  the  momentum  thickness  Reynolds  number  at  the  end  of  transition,  ReeT, 
by 

Re/r  =  75  Re|r  .  (2.2) 

Thus,  we  have  three  Reynolds  numbers:  Re  enters  the  asymptotic  develop¬ 
ment  whereas  processing  of  experimental  findings  is  based  on  Re^  and  Re^- 
Some  comments  are  due  at  this  point.  Firstly,  for  a  flat  plate,  the  momentum 
thickness  Reynolds  number,  Re#,  can  be  expressed  through  the  displacement 
thickness  Reynolds  number,  Rr-y ,  as 


Re#  =  0.3859  Re^  .  (2.3) 

An  analogous  relationship  holds  for  any  airfoil  section,  with  the  numerical 
factor  on  the  right-hand  side  depending  on  its  shape.  Secondly,  within  the 
framework  of  the  high  Reynolds  number  approach  the  transition  length  Ip  is 
assumed  to  be  much  shorter  than  the  characteristic  distance  L*.  Asymptot¬ 
ically,  as  e  tends  to  zero,  (2.2)  becomes 

Re/T  =  22.81  Ref  (2.4) 

with  the  coefficient  determined  by  using  (2.3).  Insofar  as 

Re<5  =  1.7208  Re  5 

for  the  Blasius  boundary  layer  on  a  flat  plate,  we  finally  have  an  expression 

Re/r  -  44.96  Rei  (2.5) 

defining  the  transition  length  Ip  through  the  characteristic  distance  L* .  An 
analogous  dependence  is  typical  of  an  arbitrary  blade  or  wing  section.  As¬ 
ymptotically,  (2.4)  and  (2.5)  are  tantamount  to  (2.2)  and  can  be  exploited  as 
empirical  correlations,  however  the  results  from  wind-tunnel  measurements 
collapse  without  appreciable  scatter  onto  a  straight  line  drawn  by  Mayle 
(1991)  from  (2.2).  For  this  reason  the  momentum  thickness  Reynolds  num¬ 
ber  seems  to  be  preferable  when  processing  experimental  data. 


13 


Let  //*  and  T*  be  the  viscosity  and  temperature,  respectively.  With  the 
subscript  oo  labelling  the  conditions  at  the  upper  edge  of  the  boundary  layer, 
the  Chapman-Rubesin  linear  law  frequently  used  in  theory  reads 


Ml 


C  =  const. 


The  triple-deck  time  ,  **  ^ ,  where  is  the  velocity  of  the  oncoming  stream, 

and  spatial  coordinates  are  scaled  in  terms  of  Reynolds  number 

through  the  small  parameter  e  introduced  in  (2.1)  and  include  the  ratio 
of  the  wall  temperature  T*  to  the  external  temperature  as  well  as 
t,£e  Mach  number  and  the  nondimensional  skin  friction  A  =  0.3321,  in 
the  incompressible  Blasius  flow.  For  an  arbitrary  boundary  layer,  the  skin 

1  /  'jp*  \  1 

friction  tw  can  be  expressed  through  A  by  means  of  rw  =  A (  jf-J  • 
The  above  scaling  defines  a  characteristic  Strouhal  number  for  free  interac¬ 
tion.  According  to  Stewartson  (1974),  Smith  (1982),  and  Kluwick  (1998), 
the  independent  variables  in  the  viscous  near-wall  sublayer  of  the  subsonic 
boundary  layer  are 


e2Cl<  A-l  (Af£,  -  1)  * 

£3Ci A“§  (M^-l)_i 
e5ci A“t  (jl^-lp 


Related  to  these  normalized  coordinates,  the  velocity  field 


(2.6a) 


(2.6b) 


(2.6c) 


u 

U* 


U* 


£3CiAi  {Ml,-  1)* 


and  the  pressure  variations 


(2.7a) 


(2.7b) 
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(2.8) 


’L-l~=e>Chi(Ml-l)-ip 

r  00^00 

obey  the  system  of  Prandtl  equations 


du  du  du 

di  +  uai  +  %, 


du  dv 
dx  +  dy  ’ 

(2.9a) 

dp  d2u 
dx  dy 2 

(2.9b) 

for  an  incompressible  fluid,  independent  of  the  Mach  number.  The  self- 
induced  pressure  p  =  p(t,x)  comes  from  viscous/inviscid  interaction  rather 
than  being  given  beforehand  by  a  potential  flow  solution.  In  the  absence 
of  perturbing  agencies,  the  pressure  rises  and  drops  are  connected  with  the 
instantaneous  displacement  thickness  —A(t,x)  by  the  interaction  law 


The  limit  condition 


dA(t,0 

dt 


(2.10) 


u  —  y—>A(t,x )  as  y  — >  oo  (2.11) 

at  the  upper  reaches  of  the  viscous  near- wall  sublayer  and  the  no-slip  condi¬ 
tion 


u  =  v  =  0  (2.12) 

at  the  body  surface  make  the  mathematical  formulation  complete. 

In  order  to  extend  the  above  consideration,  we  may  drop  the  Chapman- 
Rubesin  viscosity  law  and  introduce  the  independent  variables  through 


-  £2rl{M\ 


2 

oo 


t. 


(2.13a) 

(2.13b) 

(2.13c) 
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instead  of  (2.6a-2.6c,).  The  nondimensional  velocity  field  and  the  pressure 
variations 


u 

U* 

v* 

U) 

P*  "Poo 

(frJJ. 


oo 

* 


2 

OO 


£Tw  {Ml  - 
M  {Ml 
e2A  {Ml 


(2.14a) 

(2.14b) 

(2.14c) 


come  in  place  of  (2.7a,  2.7b)  and  (2.8),  respectively.  Here  rw  and  the  ratio 
j*2-  depend  on  a  particular  boundary  layer  in  question.  The  scaling  in  terms 
of £  as  well  as  M^-dependence  remain  intact,  however  the  normalization  that 
accounts  for  the  temperature  factor  varies  with  the  type  of  the  boundary  layer 
under  consideration.  Nonlinear  stability  of  boundary  layers  for  disturbances 
of  various  sizes  was  discussed  by  Smith  (1979a). 


The  length  l*  of  the  nonlinear  interaction  process  governed  by  the  system 
of  Prandtl  equations  can  be  estimated  as  l*  ~  x*  with  x*  defined  in  equations 
(2.6b)  and  (2.13b).  Hence  it  follows  that 

l*  ~  e3L*  =  Re"  I  L* 


in  view  of  (2.1).  The  Reynolds  number  Re,,  calculated  with  this  characteristic 
length  is  expressed  by  virtue  of 

l*  5 

R.e,j  =  Re  ~  Re§  (2.15) 

Ju 

in  terms  of  the  Reynolds  number.  Thus,  the  transition  length  Ij,  in  (2.4) 
exhibits  the  same  scaling  as  the  interaction  length  !*  in  (2.15).  In  fact, 
can  be  identified  with  l*  because  evidently  <  l*.  The  result  obtained 
goes  far  beyond  the  scope  of  the  triple-deck  theory  which  has  been  originally 
invented  as  applied  to  laminar  boundary  layers.  The  identity  of  the  two 
lengths  implies  that  the  asymptotic  approach  extends  to  cover  the  regime  of 
laminar /turbulent  transition  (except  when  the  intermittency  becomes  close 
to  1). 
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Simple  physical  reasoning  provides  arguments  substantiating  this  infer¬ 
ence.  The  triple-deck  scaling  of  time  and  spatial  coordinates  is  typical  of  vis- 
cous/inviscid  interaction.  On  the  other  hand,  the  same  type  of  interaction  is 
known  to  control  the  boundary-layer  eigenmodes  in  the  Tollmien-Schlichting 
spectral  range  centered  about  the  lower  branch  of  the  neutral  stability  curve 
(Smith  1979b;  Zhuk  &  Ryzhov  1980).  Thus,  the  scaling  introduced  above 
spesifies  the  period  and  wavelength  of  oscillation  cycles  within  the  Tollmien- 
Schlichting  wave  packets.  The  zone  of  transitional  flow  was  recorded  in  nu¬ 
merous  wind-tunnel  tests  to  be  filled  with  vortex  spots  (wave  packets)  that 
amplify  exponentially  fast  and  turn  into  turbulent  spots  coalescing  with  each 
other  and  thereby  sharply  increasing  intermittency  at  later  stages  of  transi¬ 
tion.  Between  these  spots  the  flow  field  remains  laminar.  So  long  as  there 
is  enough  space  between  neighbouring  spots  they  develop  as  if  they  were 
travelling  in  laminar  surroundings.  This  is  the  reason  why  the  transition  and 
free  interaction  lengths  have  equal  scaling  in  terms  of  the  Reynolds  num¬ 
ber.  Vigorous  growth  of  the  pulsation  amplitudes  within  the  spots  does  not 
change  the  scaling  in  (2.15)  but  can  trigger  erratic  short  wavelength  wiggles 
to  destroy  the  original  well-organized  vortex  structures. 


3  Numerical  evidence 

The  last  statement  needs  to  be  reinforced  computationally.  To  this  end,  let 
us  consider  two  mechanisms  of  the  vortex-spot  production. 

Wave-packet  excitation  by  a  convected  vortex.  The  first  mecha¬ 
nism  comes  about  by  virtue  of  the  potential  vortex/surface  roughness  inter¬ 
action.  Strong  similarities  in  the  transition  process  between  wake-induced 
turbulent  patches  on  turbomachinery  blades  and  artificially  triggered  turbu¬ 
lent  spots  in  the  flat  plate  boundary  layer  underlie  this  theoretical  model.  As 
mentioned,  these  similarities  were  exploited  by  Gostelow,  Hodson  &  Walker 
(1999)  and  Gostelow  &  Thomas  (2005)  in  their  wind-tunnel  tests  to  simulate 
the  transition  phenomena.  Certainly,  natural  transition  on  turbomachinery 
blades  is  three-dimensional.  However,  in  Gostelow  &  Thomas  (2005)  wakes 
were  introduced  into  the  flow  by  a  tapered  rod,  of  average  radius  4.5  mm, 
contilevered  from  a  disk  having  axis  0.5  m  upstream  of  the  leading  edge  of 
the  thin  plate.  The  rod  was  mounted  at  a  radius  of  270  mm  on  the  disk 
and  was  rotated  at  an  angular  velocity  of  60  rpm,  resulting  in  the  introduc¬ 
tion  of  two  dissimilar  wakes  per  second.  Centerline  phase-averaged  velocity 
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Figure  1:  Schematic  of  vortex/hump  interaction. 


traces  were  obtained  starting  from  the  aforementioned  similarities  between 
the  wake- induced  turbulent  patches  and  artificially  excited  turbulent  spots. 
In  fact,  the  covected- vortex  model  relates  to  conditions  typical  of  the  experi¬ 
mental  set-ups  rather  than  to  the  general  case  for  which  no  wind-tunnel  data 
are  available. 

Figure  1  gives  an  idea  of  the  process  and  makes  clear  the  triple-deck  pat¬ 
tern  of  the  velocity  field  with  pertinent  scalings  involved.  We  point  out  briefly 
the  most  important  features  associated  with  the  convected  vortex.  Asymp¬ 
totically,  as  £  — >  0,  the  flow  in  the  near-wall  sublayer  obeys  the  system  of 
Prandtl  equations  (2.9a,  2.9b).  However  an  additional  term  is  included  in 
the  expression  (2.10)  for  the  excess  pressure  to  take  into  account  the  vortex- 
induced  contribution.  The  limit  condition  (2.11)  with  the  instantaneous  dis¬ 
placement  thickness  entering  the  right-hand  side  also  contains  a  correction 
pv  (t)  for  the  vortex-induced  pressure.  This  inhomogeneous  limit  condition 
at  the  upper  edge  of  the  near-wall  sublayer  should  be  supplemented  with  the 
inhomogeneous  limit  condition  at  the  entry,  as  x  — >  — oo,  to  the  sublayer. 
The  derivation  of  the  inflow  condition  rests  upon  a  detailed  analysis  of  the 
vortex  motion  in  a  region  upstream  of  the  local  roughness.  Matching  with 
the  Stokes  sublayer  located  close  to  the  wall  in  this  region  shows  that  the 
instantaneous  displacement  thickness  tends  to  zero  as  x  — >  —  oo,  while  the 


vortex-induced  velocity 


u—* 


y 

y1 

/  ,  .3  eXP 

/-oo  (t  —  t)  2 

.  4(*-r). 

Pv  (t)  d,T 


(3-1) 


provides  the  time-dependent  contribution.  In  consequence, 


u-y  ^pv(t)  +  U  (£,  y) 

as  x  — ►  — oo.  The  no-slip  condition  (2.12)  applies  at  the  body  surface  with  a 
small  hump  or  dent 


V  —  Vw  —  ag  (x)  (3.2) 

centered  around  the  origin. 

Following  Burggraf  &  Duck  (1982)  and  Duck  (1987)  the  boundary-value 
problem  posed  is  recast  by  introducing  the  shear  stress  r(t,x,y)  as  a  new 
desired  function.  Computing  the  Fourier  spectral  distribution 

r  (t,  k,  y)  =  dU  V ^  =  —  J  t  (t,  x,  y )  e~lkxdx 

and  then  inverting  it  to  the  physical  space 

r  (t,  x,  y)  =  V ^  ^  J  t  (t,  k,  y)  elkxdk 

play  a  dominant  role  in  an  iteration  procedure  devised  for  an  equation 

dr  dr  dr  d2r 

dt  dx  dy  dy 2 

that  arises  upon  differentiating  (2.9b)  with  respect  to  y.  Our  main  concern  is 
with  the  onset  of  erratic  wiggles  in  tw  —  r  (t,  x,  0)  accompanied  with  filling  up 
the  spectrum  tw  =  t  (t,  k,  0)  with  t  increasing.  A  code  with  FFT  algorithm 
has  been  tested  in  Ryzhov  &  Savenkov  (1989)  and  Ryzhov  &  Timofeev  (1995) 
when  analyzing  the  short-scaled  large-amplitude  distortions  appearing  in  a 
few  central  soliton-like  cycles  of  the  wave  packets  which  can  be  identified 
with  the  vortex  spots  in  laminar /turbulent  transition.  The  geometry 

g  =  ne-]kl  (3.3) 

1  +  x* 
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assumed  in  the  computation  specifies  a  hump  mounted  on  an  otherwise  flat 
plate,  accordingly,  the  height  of  the  obstacle  in  (3.2)  being  a  =  1. 

For  better  visualization  of  oscillation  properties  a  quantity 


7"  zzr  7"  —  1  — 

1  w  1  w 


m^y) 

dy 


(3.4) 


is  meant  by  r)'  in  what  follows.  This  quantity  equals  r)  =  tw  ••  1  for 
the  steady  flow  set  in  prior  to  the  vortex  motion  that  induces  in  line  with 
(3.1)  an  additional  component  U  ( t,y )  of  velocity.  The  hump  fixed  by  (3.3) 
introduces  moderate  disturbances  in  the  velocity  field  of  the  oncoming  steady 
stream  without  separation  accompanied  by  a  bubble  with  reversed  flow  (for 
separation  to  occur  the  hight  of  the  hump  should  be  larger).  The  resulting 
boundary-value  problem  provides  the  simplest  model  of  gas-turbine-engine 
environments  where  vortical  disturbances  interact  with  small  imperfections  of 
the  turbine-blade  surface  (especially  covered  by  the  thermal  barrier  coating) . 
Recent  wind-tunnel  tests  by  Pinson  &  Wang  (2000)  lend  credence  to  the 
applicability  of  the  transition- length  correlation  to  rough  surfaces. 

Figures  2a,  b,  c  show  the  origin  and  rapid  development  of  a  wave  packet 
induced  by  a  convected  vortex  interacting  with  the  hump  shaped  by  (3.2), 
(3.3).  A  value  p™  =  max  (|p„|)  =  3  determines  a  fairly  strong  intensity  of 
the  vortex.  At  t  =  0,  the  excess  shear  stress  t%  is  associated  predominantly 
with  the  steady  perturbations  from  the  hump,  the  formation  of  unsteady 
oscillations  downstream  is  barely  discernible.  The  explosive  character  of  the 
highly  modulated  wave  packet  propagation  becomes  clear  at  t  =  1.5  when  a 
narrow  soliton-like  cycle  appears  in  the  central  part  of  the  signal.  Figure  2c 
illustrates  the  onset  at  t  =  2.5  of  short-scaled  erratic  wiggles  that  distort  a 
portion  of  two  central  cycles  of  the  wave  packet  with  positive  rdw.  The  soliton- 
like  portions  of  the  same  cycles  with  negative  remain  almost  intact  at  this 
spiky  stage  typical  of  transition. 

Figures  3a,  b,  c  display  the  spectral  content  of  disturbances  shown  in 
Figures  2a,  b,  c.  The  first  of  them  presents  the  wavenumber  spectrum  of 
almost  steady  perturbations  triggered  by  the  hump,  with  the  peak  in  |r^| 
located  approximately  at  k  —  ks  =  1.19.  The  second  plot  exhibits  the 
spectrum  featuring  the  wave  packet  at  t  =  1.5.  The  local  maximum  at  k  =  ks 
becomes  much  smaller  compared  to  the  local  peak  at  about  k  =  km  —  2.75 
that  corresponds  to  the  wavenumber  of  the  fastest  growing  linear  eigenmode. 
The  second  local  maximum  at  k  «  2 km  =  5.5  forks  into  two  small  mounds 
divided  by  a  shallow  local  minimum  in  Ir^l.  The  third  local  smoothed-out 
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maximum  arises  at  k  «  3 km  =  8.25.  One  more  local  maximum  is  scarcely 
visible  close  to  k  «  Akm  =  11.0.  Thus,  the  filling-up  of  the  wavenumber 
spectrum  basically  follows  a  weakly  nonlinear  scenario  in  spite  of  the  fact 
that  the  signal  in  Figure  2b  is  essentially  nonlinear.  To  the  contrary,  the 
filling- up  of  the  distant  parts  of  the  spectrum  in  Figure  3c  is  at  the  heart  of 
the  short-scaled  wiggles  corrupting  two  central  cycles  of  the  signal  in  Figure 
2c.  Here  the  first  two  peaks  in  |  |  become  almost  equal  whereas  the  other 

local  maxima  attain  the  magnitudes  comparable  to  these  peaks.  The  signal  as 
a  whole  nears  to  transitioning  to  a  turbulent  spot.  This  process  embracing  all 
the  vortex  spots  in  the  boundary  layer  eventually  brings  up  the  intermittency 
close  to  1. 

Acoustic  wave/surface  roughness  interaction.  The  second  mecha¬ 
nism  of  the  vortex  spot  production  comes  about  by  virtue  of  sound  scattering 
into  Tollmien-Schlichting  eigenmodes  by  small  streamwise  variations  in  sur¬ 
face  geometry.  This  mechanism  is  even  more  important  as  applied  to  aircraft 
wings  rather  than  turbomachinery  blades.  Time-harmonic  oscillations  were 
studied  by  Goldstein  (1985)  and  Goldstein  &  Hultgren  (1989).  The  closest 
resemblance  to  the  wave-packet  excitation  by  a  convected  vortex  shows  the 
propagation  of  a  sound  pulse  impinging  against  a  local  obstacle  located  on  an 
otherwise  smooth  surface  (Ryzhov  &  Timofeev  1995).  However  there  exists 
a  difference  between  the  two  problems  because  a  fluid  in  the  potential  flow 
domain  cannot  be  regarded  anymore  to  be  incompressible  even  if  velocities 
are  small.  The  Mach  number  and  temperature  ratio  dependencies  should 
be  incorporated  into  a  set  of  similarity  parameters  in  a  manner  specified  in 
(2.6a-2.6c),  (2.7a,  2.7b)  and  (2.8),  or  (2.13a-2.13c)  and  (2.14a-2.14c).  The 
system  of  incompressible  Prandtl  equations  (2.9a,  2.9b)  still  controls  the  flow 
field  in  the  viscous  near-wall  sublayer  where  the  excess  pressure  related  to  the 
instantaneous  displacement  thickness  obeys  (2.10)  with  an  additional  term 
accounting  for  sound  scattering  at  subsonic  velocities.  What  is  more,  the  form 
/  (t)  of  the  sound  pulse  can  be  defined  in  such  a  way  as  to  render  it  identical 
to  a  function  pv  ( t )  that  was  used  above  to  prescribe  the  vortex-induced  pres¬ 
sure.  This  analogy  enables  the  wave-packet  excitation  by  a  convected  vortex 
to  be  experimentally  simulated  by  means  of  acoustic  generators. 

Instead  of  applying  this  analogy  we  put  in  the  computation 


/  =  0  for  t  <  0 
/  =  O.lart  exp  (1  —  tot)  for  t  >  0 
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